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Abstract 

We develop an algebraic approach to studying the spectral proper- 
ties of the stationary Schrodinger equation in one dimension based on 
its high order conditional symmetries. This approach makes it possible 
to obtain in explicit form representations of the Schrodinger operator 
bynxn matrices for any n € N and, thus, to reduce a spectral prob- 
lem to a purely algebraic one of finding eigenvalues of constant n x n 
matrices. The connection to so called quasi exactly solvable models is 
discussed. It is established, in particular, that the case, when condi- 
tional symmetries reduce to high order Lie symmetries, corresponds to 
exactly solvable Schrodinger equations. A symmetry classification of 
Schodinger equation admitting non-trivial high order Lie symmetries 
is carried out, which yields a hierarchy of exactly solvable Schrodinger 
equations. Exact solutions of these are constructed in explicit form. 
Possible applications of the technique developed to multi-dimensional 
linear and one-dimensional nonlinear Schrodinger equations is briefly 
discussed. 
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I. Introduction 



Basic motivation for introducing conditional symmetries (the term 'condi- 
tional symmetry' was suggested for the first time by Fushchych 0-0) 
was a necessity to find a symmetry background of a quickly growing vari- 
ety of exact solutions of nonlinear partial differential equations that could 
not be obtained within the framework of the classical Lie approach. An 
intensive search of such solutions was begun independently and almost si- 
multaneously by Fushchych with collaborators (see, ]|, Q and references 
therein), Clarkson & Kruskal ('the direct reduction method' f|), Olver 
& Rosen au ('non-classical reduction' 0) and Winternitz & Levi ||. A 
number of examples of nonlinear partial differential equations in two, three 
and even four dimensions having non-trivial conditional symmetries is grow- 
ing rapidly. In particular, it has been established by Fushchych, Zhdanov 
and Revenko Pi— jlql that such fundamental equations of the modern quan- 



tum field theory as the four- dimensional nonlinear d'Alembert, Dirac, Levi- 
Leblond, Maxwell and Yang-Mills equations possess infinite conditional sym- 
metries, while their Lie symmetries are finite only. 

On the other hand, much less attention is devoted to the study of con- 
ditional symmetries of linear differential equations (though the first example 
of conditional symmetry has been obtained by Bluman and Cole for the 
one-dimensional linear heat equation fljj]). In view of the role played by 
conditional symmetries in the theory of nonlinear differential equations one 
can expect that application of these to linear equations will also be rich in 
results. In the present paper we establish the rather unexpected (at least 
for the author) fact that conditional symmetries can be effectively applied to 
study spectral properties of the stationary Schrodinger equation 

i> xx = {e + V{x))ij. (1) 

In particular, we will prove that it is conditional symmetry that is respon- 
sible for a phenomenon of so called 'quasi exact solvability' of some specific 
class of equations ([I]) [|T^]-[|I7|. And what is more, the case when conditional 
symmetries are equivalent to Lie symmetries will be shown to yield exactly 
solvable Schrodinger equations. 

All principal approaches to a construction of (quasi) exactly solvable mod- 
els (apart from the specific ways of an implementation of these) are based 
on a possibility to construct a basis wi(x), U2(x), . . .,u n (x) of some invariant 
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space V n of the Schrodinger operator S = d% — V(x). This means that there 
should exist constant n x n matrix ||Ajfc|| such that the following conditions 
are fulfilled: 



Suj (x) 



3 = 1,2, 



n. 



(2) 



fc=i 



Given such functions procedure for calculating the spectrum of 

the Schrodinger operators (or, more precisely, a part of the spectrum) is 



completely algebraic. Let a 3 



complete system of eigenvectors of the 
Ai, . . . , A m be their eigenvalues, namely 



'n x n) 



3=},... 
matrix A 



, m, m < n be a 



I A 



■jk\\j,k=\ 



n 

E 



K jk a] = X t a l k , j 



1,2, 



, m. 



Then, the function 



x 



and 



(3) 



(4) 



3=1 



is easily seen to satisfy the equation (P with e = \k under arbitrary k = 
1,2,..., m. 

Saying it another way, after being restricted to a linear space V n with basis 
functions u\(x), u 2 (x), . . -,u n (x) the Schrodinger operator becomes a matrix 
operator. Thus, a reduction of a differential operator to a matrix operator 
takes place. But such a procedure is quite a common routine in the theory of 
Lie symmetries of differential equations. Indeed, if we restrict a partial dif- 
ferential equation having N independent variables to a subset of its solutions 
invariant under a one-parameter subgroup of the Lie group admitted by the 
equation in question, then it is reduced to a partial differential equation with 
N — 1 independent variables. Such a procedure is called symmetry reduction 
of differential equations (for more detail, see e.g. 0, [31| |19|). Taking N = 1 



(the case of ordinary differential equation) we obtain as a reduced equation 
a differential equation with N = 0, i.e. an algebraic equation! 

One of the main aims of the present paper is to show that the idea of 
symmetry reduction, when formulated in an appropriate way, can be applied 
effectively to an algebraization of the problem of describing spectrum of the 
Schrodinger operator. 
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As mentioned above classical Lie symmetries of partial differential equa- 
tion do not give all possible reductions. More general symmetries responsible 
for a possibility of reducing the order of differential equations are conditional 
symmetries. Roughly speaking, the necessary and sufficient condition provid- 
ing a possibility to reduce a number of variables in a given partial differential 
equation is a requirement of conditional invariance ( ||20|| ). It will be estab- 
lished that a similar situation takes place for the Schrodinger equation ([!]). 
Symmetries providing reducibility of differential equation ([[]) to a system of 
algebraic equations of the form (0) are exactly the high order conditional 



symmetries introduced independently by Zhdanov & Fushchych [^T], ^2 
and Fokas & Liu || (see also H-||). 



It should be emphasized that considerations of the present paper are 
purely algebraic. The method of conditional symmetries making it possible 
to study spectral properties of the Schrodinger operator S gives no informa- 
tion about analytical properties of the corresponding eigenf unctions. In each 
specific case such properties of the eigenfunctions obtained as square inte- 
grability, asymptotic behavior, singularities etc. should be studied indepen- 



dently (see, e.g. |p7[). The reason is that this method (and group-theoretical, 
symmetry methods in general) exploits algebraic properties of the solution set 
of equation (D) (or its part) as a whole and, roughly speaking, is independent 
of analytical properties of specific solutions. 



II. Conditional symmetry of the Schrodinger 
equation 

Consider the n-th order differential operator 

n 

Q = J2 qj (x)dl, (5) 
j=o 

where <9° = 1, d x +1 = 4- and functions qj(x) are supposed to be indepen- 
dent of e. 

Following [28, 29] we say that equation ([I]) is conditionally invariant with 



respect to the operator Q if the following operator identity holds: 

[Q, d 2 x -(e + V{x))} = RQ + P (dl-(e + V{x)). (6) 
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Here [Qi, Q2] = Q1Q2 — Q2Q1, R, P are some first and n-th order 
differential operators, correspondingly. The above operator equality should 
be understood in the following way: the differential operators in the left- and 
right-hand sides give the same result when acting on arbitrary (n + 2)-times 
continuously different iable function f(x). 

Provided R vanishes, condition (|6]) is nothing else but a criterion for 
equation (P to be invariant with respect to the operator Q. In such a case, 
the operator Q is a generalized (high order) Lie symmetry operator. But 
given a condition R 7^ 0, the operator Q corresponds to high order conditional 
symmetry of the Schrodinger equation ([!]). 

It is easy to see that if an equation is conditionally invariant with respect 
to the operator Q, then it is conditionally invariant with respect to the opera- 
tor q(x)Q with an arbitrary sufficiently smooth function q(x). Consequently, 
without loss of generality we can suppose that in q n (x) = 1 and consider 
differential operators of the form 

n-l 

Q = 3? + E (7) 
3=0 

As the coefficients of the operator Q do not depend on e, equality ([]) is 
only possible if R = r(x) and P = with some sufficiently smooth function 
r(x). Consequently, the condition (^) is rewritten to become 

[Q, d 2 x - V{x)\ = r(x)Q. (8) 

We call the Schrodinger operator reducible if there exist linearly-indepen- 
dent functions Ui(x), . . . , u n (x) and constants Aj k such that the conditions (g) 
are fulfilled. Let us note that this terminology is justified both from the point 
of view of the classical representation theory and of the symmetry analysis of 
differential equations. Indeed, conditions (0) mean that the representation 
space of the operator S contains an invariant subspace and, consequently 
the representation is reducible. On the other hand, conditions (0) ensure the 
reduction of the differential equation flU) to a system of algebraic equations. 
We will prove an assertion which shows that this is not a simple coincidence 
but a fundamental fact having a natural symmetry interpretation. 

Theorem 1 The Schrodinger operator S = d 2 x — V(x) is reducible if and 
only if there exists an n-th order differential operator Q of the form ^) such 
that equation is conditionally invariant with respect to Q. 
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Proof. The necessity. Let the operator S = d\ — V(x) be reducible. Then, 
the conditions @ hold. As functions Uj(x) are linearly independent, they 
form a fundamental system of solutions of some n-th order linear ordinary 



differential equation [p0[ . 



We recall that fundamental system of solutions of an ordinary differential 
equation is a maximal set of its particular solutions such that any smooth 
solution can be represented as a linear combination of these. Provided the 
order of the ordinary differential equation in question is equal to n, any n lin- 
early independent solutions of it form a fundamental system. Furthermore, 
having a fundamental system of solutions we can reconstruct the correspond- 
ing ordinary differential equation within a multiplication by a function r(x). 
Consequently, if we fix the coefficient of the n-th order derivative to be equal 
to 1, then this equation is unique. 

Thus, there exists the n-th order differential equation 

n-l 

u [n \x) + qj(x)u U) (x) = 0, (9) 

3=0 

such that the functions Uj(x) form a fundamental system of its solutions. 

We will prove that equation (|l|) is conditionally invariant with respect to 
the operator 

n-l 
3=0 

By force of relations @ the following equalities hold: 

[Q, d 2 x - V(x)] Uj (x) = Q{{dl - V(x))u 3 (x)} 

-{dl - V(x)){Q Uj (x)} = Q \J2 Ai*«*( 

U=i 

for any j = 1, 2, . . . , n. 

Thus, the functions Uj(x) satisfy an ordinary differential equation 

[Q, d 2 x -(e + V(x))}u(x) = 0, (10) 

whose order is easily established to be equal to n. Consequently, its fun- 
damental system of solutions consists of n functions. Hence, we conclude 
that the functions Uj(x) form a fundamental system of solutions of (|T0"D. As 
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an ordinary differential equation is determined by its fundamental system 
uniquely within a multiplication by a function r(x), the relation hold 

[Q, d 2 x -(e + V(x))]=r(x)Q, 

which is the same as what was to be proved. 

The sufficiency. Let the Schrodinger equation ([]]) be conditionally invari- 
ant with respect to the operator (|7|), which means that the condition (H) is 
fulfilled. Consider an equation 

Qu(x) = (d: + u{x) = (11) 

as an ordinary differential equation for a function u(x). Clearly, its general 
solution is represented in the form 

n 

u [x) = CjUjix), (12) 

3=1 

where Cj are arbitrary constants and u±(x), . . . , u n {x) is a fundamental sys- 
tem of solutions of n-th order ordinary differential equation ([II]) . 

From the condition (§) it follows that the Schrodinger operator S = 
d 2 — V(x) is a symmetry operator for the equation (|ll|), i-e. it transforms 
each solution of equation ( ITT] ) into another solution of the same equation. 



Consequently, for any j = 1,2, ... ,n the function Uj(x) = (d^ — V(x))uj(x) 
satisfy (|TTD. But by definition the fundamental system of solutions of n-th 
order ordinary differential equation forms a maximal set of its linearly inde- 
pendent solutions, which means that any solution can be represented as a 
linear combination of functions Uj(x). Thus, there exist such constants Ajk 
that functions Uj(x) satisfy relations (0), whence it follows that the corre- 
sponding Schrodinger equation is reducible. The theorem is proved. 

Note that the proof of theorem is, in fact, independent of the specific form 
of the Schrodinger operator S = d 2 x — V(x). It is straightforward to generalize 
Theorem 1 to the case of an arbitrary iV-th order differential operator 

N 

S = Y,fi(?)8L. (is) 

We give the corresponding assertion without proof. 
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Theorem 2 The operator S in flT3i) is reducible if and only if there exists n- 
th order differential operator Q of the form (Qj such that equation Sip{x) = 
is conditionally invariant with respect to Q. 

To illustrate the above statement we consider two examples. 
Example 1. Consider the harmonic oscillator Schrodinger equation 

^ xx = {e + x 2 )^. (14) 

As a direct check shows the n-th order differential operator 

Q = (d x -x) n (15) 

satisfies the following commutation relation: 

[Q, d 2 x -{e + x 2 )} = 2nQ 

(the easiest way to prove the above formula is to use the mathematical in- 
duction method). 

Consequently, equation (|i"4]) is conditionally invariant with respect to the 
operator Q and we can apply Theorem 1. Integrating equation Qtp(x) = 
yields a basis of the invariant space V n of the Schrodinger operator d 2 — x 2 

2 2 2 2 
x x q x i x 

P 2 TP 2 TP 2 X P 2 

It is readily seen that the above functions satisfy relations (0) with V(x) = 
x 2 . Calculating eigenvalues (Xj) and eigenvectors (a- 7 ) of the corresponding 
matrix ||A^|| we obtain exact solutions of the Schrodinger equation ( |i~4"D with 
e — Ai, . . . , X m in the form (^). 

Example 2. Let us generalize the previous example as follows. We are looking 
for the Schrodinger equations ([!]) conditionally invariant with respect to the 
n-th order operator which can be represented as a power of the first order 
differential operator, i.e. 

Q=(a(x)d x + b(x)) n+1 . (16) 

By an appropriate transformation of the dependent and independent vari- 
ables 

z = F(x), <p(z) = e~I G{x)dx ^P(x), 
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we can transform the operator Q as follows: 



Q = d[ 



tn+l 
z 



After rewriting the initial Schrodinger equation in the new variables 

z, ip(z) we get 

f(z)(f zz + g(z)<p z + (h(z) - e)(p = 0, 

where 

f(z) = (F'(x)) 2 , g{z) = F"(x) + 2F'(x)G(x), 
h(z) = -V(x) +G'(x) + G 2 (x). 

Commutation relations (f|) now read as 

[8p\ f(z)d 2 z + g(z)d z + h(z) -e}= r(z)d: +1 . (17) 

Computing the commutator in the left-hand side (which is a simple ex- 
ercise in differential calculus) and equating coefficients of the linearly inde- 
pendent operators d{ we conclude that the equation ([H]) is consistent if and 
only if the functions f,g,h are polynomials in z of the following form: 

h(z) = A -n(B 2 + (n-l)C 3 )z + C 4 (n~l)nz 2 } 
g{z) = B + B x z + B 2 z 2 + 2C 4 (1 - n) z 3 , 
f{z) = C + C 1 z + C 2 z 2 + C 3 z 3 + C A z\ 

where A , B , B\, . . . , C4 are arbitrary constants. 

Returning back to the initial variables x, u(x) we get the necessary and 
sufficient conditions for the Schrodinger equation (H) to be conditionally in- 
variant with respect to an operator belonging to the class ( |TED 

-V + G' + G 2 = A - n (B 2 + (n - 1)C 3 ) F + C 4 (n - 1) nF 2 , 
F" + 2F'G = B + BiF + B 2 F 2 + 2C A (1 - n) F 3 , 
(F') 2 = C + CtF + C 2 F 2 + C 3 F 3 + C 4 F\ 

whence we derive the form of the potential V(x) 

y, v ^0 + V\U + v 2 uj 2 + V 3 UJ 3 + v^u* 
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v = AB 2 - 16A C + 8SiC - 8B Ct + 3C ± 2 - 8C C 2 , 
vt = SB Bi + 16B 2 Co(n+l) -16A C 1 -16B C 2 + ^C 1 C 2 

+ !6B 2 C n + 8C C 3 (2n 2 - 2n - 3), 
v 2 = 45 1 2 + 85 5 2 + 85 2 C 1 (2n + l) - 16^2-85^2 + 4C 2 2 

-2AB C 3 + 2CiC 3 (8n 2 - 8n - 3) - 16C C 4 n(n + 2), 
v 3 = 8B t B 2 - 16A C 3 - I6S1C3 - 16£(A(7i + 1) + 16B 2 C 2 n 

-16B C 4 n - 8CiC 4 (2n 2 + 2n - 1) + AC 2 C 3 (An 2 -An + 1), 
w 4 = 4 J B2 2 + 8 J B2C3(2n-l)-16A C4-8 J B 1 C 4 (2n + l) + 

+C 3 2 (16n 2 - 16n + 3) + 8C 2 C 4 (1 - 2n 2 ) 

and of the function 

25 - Ci + 2(Bi - C 2 )o; + (25 2 - 3C 3 )^ 2 - 4C 4 no/ 



G(x) 



In the above formulae u(x) is an elliptic function determined by the 
quadrature 

u>(x) 

dr 

x. 



VC + C x r + C 2 r 2 + C 3 t 3 + C 4 r 4 

Furthermore, exact solutions of the Schrodinger equation with the poten- 
tial (pl|) read as 

n 

=eS G ^ dx Y^a j uj(x)\ 

3=0 

where a = (ao, a\, . . . , a n ) is an eigenvector of some (n+lxn + 1) constant 
matrix whose entries are linear combinations of the parameters Aq,Bq,Bi,. . ., 
C 4 (we omit the corresponding formulae). 

Thus, we arrived at the nine-parameter family of quasi exactly solv- 
able Schrodinger equations obtained by Turbiner and Shifman within the 
framework of their Lie algebraic approach [T5|, |TB[ and by Ushveridze by 
means of a more general analytic approach. A detailed account of proper- 
ties of the Schrodinger equation with potentials (18) can be found in the 
monograph |I7 |. We restrict ourselves to noting that if we choose in the 
above formulae B 2 = C 3 = C 4 = 0, then the potential V(x) does not de- 
pend on n (the order of the operator Q) and, consequently, the corresponding 
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Schrodinger equation is exactly solvable. Thus, the well-known six-parameter 
family of exactly solvable Schrodinger equations is obtained. In particular, 
choosing C\ = C 2 = C 3 = C 4 = 0, C = 1 yields the harmonic oscillator 
Schrodinger equation fll4|). 

Now let us pass from particular examples to the general case in order to 
examine which constraints are imposed on the coefficients of the operator (0) 
by the requirement of conditional invariance In order to compute the 
commutator on the left-hand side of (§) we use the following identity: 

[d k x ,f(x)} = Y / C{f^\x)di, keN, (19) 

where = k\(J\(k — j)!) -1 are binomial coefficients, which is established by 
the mathematical induction method with the help of the evident identity 

[8 k x + \ f(x)] = d x [d k x , f(x)]+f(x)d k x . 
Taking into account formula flT9| ) we rewrite relation (^) as follows: 

n— 1 n— 1 i— 1 n— 1 

j=0 i=l i=0 j=0 

= r{x) U + ^q^. 

Comparing the coefficients of d™ on the left- and right-hand sides of the 
above equation we conclude that r(x) = —2q[ n _ x . Comparing the coefficients 
of the linearly independent operators d x , d x ,. . ., d x ~ l we arrive at the fol- 
lowing system of nonlinear ordinary differential equations for the functions 
q (x), qi(x),. . ., q n -i(x), V(x): 

H-i + q" + CI V<^» + J2 q t C{ V<*-» - 2 q ' n _, q 3 = 0, (20) 

i=j+i 

where j = 0, 1, . . . , n — 1 and by convention q_i = f 0, q n = f 1. 

Thus, we have n equations for n + 1 functions, which means that the 
system ( ^0[) is under-determined. As an immediate consequence of this fact 
we conclude that any Schrodinger equation ([!]) is reducible. Indeed, fixing 
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in an arbitrary way a function V = V(x) yields a second-order system of 
n ordinary differential equations for n functions qo(x), qi(x),..., q n -i(x). 
Each solution of such a system gives rise to an operator Q satisfying by 
construction condition Consequently, the conditions of Theorem 1 can 
be fulfilled with any choice of the potential V(x). 

Let us demonstrate how the results obtained can be used to study the 
spectral properties of the Schrodinger equation. Remarkably, to this end we 
do not need an explicit form of solution of the system of nonlinear ordinary 
differential equations (p0|) . It suffices to know initial values of the functions 
qj{x) and of their first derivatives q'j(x) at some point x = xq G R. We 
denote these as follows 

qj(x ) = A j+ i, q'j(x ) = B j+ i, j = 0,1,..., n-1. (21) 

All the information about spectral properties of the Schrodinger operator 
restricted to an invariant space V n with basis elements u\(x) ,u 2 (x), . . -,u n (x) 
is contained in the matrix 1 1 A ^ ^ 1 1 , which determines a transformation law (0) 
for the functions Uj(x) with respect to the action of the Schrodinger operator 
S. 

Let Q be a differential operator of the order n satisfying condition (H) 
with some choice of the function V(x). Then, by force of Theorem 1, an 
invariant space V n of the corresponding Schrodinger operator 5* is spanned 
by the fundamental system of solutions of the n-th order ordinary differential 
equation 

n-1 

u (n) (x) + Y, qj(x)u U) (x) = 0. (22) 

3=0 

We denote this system as {ui(x), u 2 {x), ...,u n (x)}. As any n linearly 
independent solutions of equation form a fundamental system of solu- 
tions, there is a freedom in the choice of the functions Uj(x). We fix these 
by imposing initial conditions. 

Let be a constant non-singular nxn matrix. Consider the following 

n Cauchy problems 

n-1 

4 n) (^) + E q,(x)v>k\x) = 0, k = 1, . . . ,n (23) 

3=0 

Uk~ ( x o) = L kj, k,j = l,...,n. 
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It is well-known from the general theory of linear differential equations 
that the above system has a unique solution, and what is more, this solution 
yields a fundamental system of solutions of equation (f22|) . 

Differentiating relations (0) n — 1 times with respect to x and excluding 
the n-th and the (n + l)-th derivatives of the functions Uj(x) with the help 
of equations we arrive at the following relations: 



k=l k=0 
n n— 1 Ti—2 

EA,rf- 2) = -E^f-ECt^^-^f, (24) 

k=l k=0 k=0 

n n—l 

EAirf - ^ = E (in-lQk -<&- q k -l ~ G^V^- 



k=l k=0 

where j = 1, . . . , n, i = 0, . . . , n — 3 and as above q_i = f 0. 



Choosing x = x in (p4|) yields 



where 



E ^-jkLki — ^ Lj^Rki, j,i — 1, 
fe=i fc=i 



= (Wa-C^V^Oro), z = l,...,n-2, 
^ fc n-i = -A k -C k n ^ n - k - l \x ), (25) 



In formulae (|25|) the index runs from 1 to n, -B/t are constants 
defined by (]21"D and is the Kronecker symbol. 

Rewriting the formulae obtained in the matrix form we have A L = L R, 
where A, L, R are constant n x n matrices with entries Ajk, Ljj., Rjki 
respectively. Hence, we derive the explicit form of the matrix A 

A = LRL~ l . (26) 

Thus we have proved the following assertion. 
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Theorem 3 Let L be an arbitrary invertible n x n matrix, and A\, . . . , A n , 
B\, . . . , B n be arbitrary constants. Then, for any choice of the function V(x) 
there exist n functions U\(x),. . ., u n (x) such that the relations (j^) hold, Kjk 
being the entries of the n x n matrix given by formulae $2$), $Zh\)- 



The above theorem has as a consequence the following important assertion 
which describes a finite part of the spectrum of the Schrodinger operator 5*. 

Theorem 4 Let Ai, . . . , X m be distinct eigenvalues of the matrix \\Rjk\\ de- 
termined by formulae (ffil). Then, for any choice of parameters Aj, Bj there 
exist linearly independent functions ipi{x), . . . , ip m {x) satisfying the Schro- 
dinger equation with e = Aj, £ = A2, ...,£= \ m , correspondingly. 

The proof follows from Theorem 2 if one takes into account that the 
matrix A is similar to R and, consequently, has the same eigenvalues Ai, A2, 
. . ., X m . The explicit form of the functions ipi(x), . . . , i[> m (x) is given by the 
formula where a? = (a{, a J 2 , ■ ■ ■ , a J n ), j — 1, . . . , m are eigenvectors of the 
matrix A corresponding to the eigenvalues Ai, A2, . . ., A m . 

Thus, using the conditional symmetry approach we were able not only to 
calculate the spectrum of the Schrodinger operator S = d\ — V[x) but also to 
construct in an explicit form a 2n-parameter family of matrix representations 
of S. 

As is seen from formulae (|25|), (p6|), there is a large freedom in choice of 
the matrix A. First, it depends on 2n arbitrary constants Aj, Bj, which fix a 



fundamental system of solutions of ordinary differential equation (|22|). Con- 
sequently, choosing specific constants Aj, Bj means fixing a representation 
space V n . Secondly, it contains an arbitrary constant nxn matrix L. The ap- 
pearance of the matrix L in the definition of A reflects a freedom in choosing 
a basis of the representation space V n . Indeed, if u\(x), . . . ,u n (x) is a basis 
of the space V n , then the functions J2k=i LjkUk(x) also form a basis with an 
arbitrary invertible constant nxn matrix ||Ljjfc||. That is why, choosing a 
specific matrix L results in fixing a basis of the representation space u\{x), 
. . ., u n (x). This freedom can be used, in particular, to obtain an orthogo- 
nal basis for V n (to this end one should apply the standard Gram-Schmidt 
orthogonalization procedure). 

But representations of the Schrodinger operator in these basises are equiv- 
alent, which is readily seen from the formula (|26|) . 
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Thus, it is established that any Schrodinger equation has n-dimensional 
invariant spaces V n with arbitrary n G N. Furthermore, we have con- 
structed the 2n-parameter family of matrix representations of the correspond- 
ing Schrodinger operator in these spaces (the formulae (p5|), fl26|)). But to 
obtain an explicit form of the basis of V n we still have 



to integrate system of nonlinear ordinary differential equations (BO), 
and 

to construct the general solution of the n-th order ordinary differential 
equation (|22|). 



We will demonstrate that using a simple trick we may avoid the necessity 
to integrate equation (|^). The said trick is based on the fact that we need 
not all solutions of (|22|) but only those which simultaneously satisfy the initial 
Schrodinger equation (|IJ). This means that we have to solve the following 
over-determined system of two ordinary differential equations: 

^xx = (e + V(x))if), 

^ n \x) + £ qAx)^\x) = 0. 

3=0 

Using the first equation and its differential consequences up to the order 
n — 2 we can exclude from the second equation all the derivatives of the func- 
tion ip of the order j > 1 and rewrite system ( p7|) in the following equivalent 
form: 

i>xx = (e + V(x))ip, 

( N N \ (28) 

E ai{x)e l d x + E &i(^)e 4 U = 0, 

\i=0 i=0 ) 

where N — [^], ctj(x), 6j(x) are linear combinations of the functions qj(x) 
with coefficients depending on V(x) and its derivatives (and, consequently, 
independent of e) and, furthermore, = 1 if n = 2iV+l and a at = 0, b]y = 1 
if n = 2N. 

The compatibility condition for the above system reads as 
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As functions dj, b{ are independent of e, coefficients of the powers of e 
should be independent of x. This requirement yields 

1) 2N + 1 ordinary differential equations 

(K a j - ka'j - hbj + Va^a^j + ^ a i a i = Ck, (30) 

i+j=k i+j=k—l 

where k = 0,1, . . . , 2N, for 2N + 2 functions ao, ■ ■ ■ , cln-i, &o, ■ ■ ■ , &iv, 
V, provided n = 2N + 1, or 

2) 2N ordinary differential equations of the form (]30f) with taking the 
values 0, 1, . . ., 2N — 1 for 2iV + 1 functions a , . . . , ajy-i, &o> • • • > &jv-i> 
V, provided n = 2N. 



In ( P0"D Cj are arbitrary constants. 

Provided conditions (|30"D are fulfilled, the compatibility condition ( p9|) 
reduces to an algebraic equation 



1) under n = 2N + 1 



2W 

£ 2iV + l + £ ^ = 0> (31) 
j=0 



2) under n = 2N 

2N-1 

„2iV 



+ £ = (32) 
i=o 

and the general solution of system ( ^8f) is given by the quadrature 

**>=M-/(ifc^H (33) 

It should be noted that equations Q3TD, ( |3"2"| ) are nothing else but charac- 
teristic equations for the matrix R defined by (|25|). Their solutions Ai, . . . , A m 
are eigenvalues and ipi(x) = ip(x)\ e =\n ■ • -, i/Jm{%) — ' l P{ x )\e=x m are eigenfunc- 
tions of the corresponding Schrodinger operator S = d 2 x — V(x). 
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III. Symmetry and exact solvability 



With the best of our knowledge the first paper, where a systematic study of 
high order Lie symmetries of the Schrodinger equations with non-vanishing 
potentials has been undertaken is the one by Nikitin, Onufriychuk and 



Fushchych In particular, for several one- dimensional exactly solvable 
models third-order symmetry operators were constructed. Furthermore, it 
was conjectured that exact integrability of the Schrodinger equation ([I]) is 
intimately connected with its symmetry properties. This conjecture has been 



confirmed in |M| , where a number of exactly solvable potentials were obtained 
by means of third-order symmetries of ([!]) and, furthermore, a method for 
integrating the corresponding Schrodinger equations was suggested. Using 
a technique developed in the previous section we will demonstrate how to 
derive exact integrability of equation (||) from its symmetry properties in the 
class of arbitrary order symmetry operators. 

The n-th order operator Q of the form (|5]) is a symmetry operator of 
equation (p]) if the condition (H) is satisfied with r(x) = 0. As coefficients 
of Q are independent of e, equality (g) is only possible when P = 0. Com- 
puting the commutator on the left-hand side of (|D with r = 0, P = and 
equating to zero the coefficients of <9™ +1 and <9™ we conclude that q n (x) = 
const, q n -i(x) = const. Consequently, any n-th order symmetry operator 
admitted by the Schrodinger equation (|l|) can be represented in the form (|^) 
with q n = C = const ^ 0, q n ~\ = Cq = const. 

First, following p4| we will consider in more detail the third-order sym- 
metry operators admitted by the Schrodinger equation (0). These are the 
lowest order symmetry operators not equivalent to usual first-order Lie sym- 
metries. The general form of a third-order symmetry operator is as follows 
(we choose C = 1) 

Q = 8l + C d 2 x + gi (x)d x + q (x), (34) 

where qo{x), qi{x) are sufficiently smooth functions to be determined from 
the invariance condition (fj) with n = 3, r = P = 0. A short computation 
yields the following expressions for the coefficients of the operator Q 

q (x) =C 2 - C V(x) - Ql (x) =C X - ^V(x), 
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where V(x) is an arbitrary solution of the third-order nonlinear ordinary 
differential equation 

-4dV + QVV -V i3) = 0. (35) 

Integrating twice the above equation we arrive at the first-order ordinary 
differential equation integrable in elliptic functions 

C 4 + 2 C 3 V(x) + 4 Ci V(x) 2 - 2 V(xf + V\xf = 0. (36) 



It was established in p4j that particular cases of almost all exactly solv- 
able potentials which can be expressed in elementary functions, such as the 
trigonometric and hyperbolic Poschel- Teller, Eckart, Kratzer potentials, po- 
tential well of finite and infinite depth, are obtained as solutions of the equa- 
tion Q3§. 

This fact imply an existence of an intimate connection between high order 
Lie symmetries and exact solvability of the Schrodinger equations. In what 
follows we will show that this is not simply a conjecture but a fundamental 
fact making it possible to classify exactly solvable models and to construct 



their exact solutions in an explicit form (see, also [24, 25]). 

It is straightforward to check that if Q is a Lie symmetry of the Schro- 
dinger equation (|l|), then Q = Q + P (d% — V(x) — e), where P is an arbitrary 
differential operator, is also its Lie symmetry. Furthermore, systems of ordi- 
nary differential equations 

{d 2 x - V{x) - e)i;(x) = 0, QVO) = 

and 

{d 2 x - V(x) - e)ip(x) = 0, Qip(x) = 

are equivalent in a sense that they have the same solutions. Making use of 
these facts we can reduce the n-th order symmetry operator Q to a first-order 
symmetry operator, the coefficients of which are N-th order polynomials in 
e (we will preserve the same designation Q for the reduced operator) 

(N \ N 

aj(x)e j \d x + ^ bj{x)e j . (37) 
j=0 J j=0 

From the invariance condition for the reduced operator 

[dl — V(x) — e, a(x, e)d x + b(x, e)} = R(x, e)(dl — V(x) — e) 
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we get a system of determining equations for the coefficients of Q 
a"{x,e) + 2b'(x,e) = 0, 

b"(x, e) + a(x, e)V\x) + 2a'(x, e) (V(x) + e) = 0, 

where primes denote differentiation with respect to x. 

Splitting the above equations by powers of e with subsequent integrating 
yields 

1 



bj(x) = —-a'j(x) + Bj, cin(x) = A 



\ 1 , ( 38 ) 



cij-i(x) = -a'j(x) — V(x)a,j(x) + - / V'(x)a,j(x) dx + Aj-\. 

In (^) A-i, Aj, Bj are arbitrary constants, j = 0, 1, . . . , JV, a_i(x) = f 0. 
Thus, the problem of describing n-th order symmetry operators of the 
Schrodinger equation is reduced to solving the recurrent relations 

1 1 r 

aj-i(x) = -a"Ax) — V(x)aAx) H — / V (x)a,j(x) dx + Aj-i (39) 

with oat(x) = An = const, a_i(x) = f 0, j = N, N — 1, . . . , 0. 

The first N relations (j = N, N — 1, . . . , 1) are solved by subsequent 
integrations yielding the expressions for the functions a (x), . . . , cln-i(x) via 
the function V(x) and its derivatives. Substituting these results into the last 
equation (j = 0) we arrive at the 2iV-th order nonlinear ordinary differential 
equation for the function V(x). It will be shown that any solution of this 
equation gives rise to an exactly solvable Schrodinger equation (|J). 

To reveal the structure of the equation in question we introduce the new 
functions Uq(x), Ui(x),. . . by the following recurrence relation: 

U j (x)=XU j ^ l (x)= (ldl-V(x) + ~d; l V'(x)} Uj-!, J = 0,1,... (40) 

where d~ l denotes integration with respect to x and U_i(x) = 1. 

Formally, a definition of functions Uj(x) contains a quadrature but inte- 
grating relations (^) successively we can get rid of it for any j = 1, 2, 3, . . .. 
Below, we adduce expressions of the functions Uj(x) for j = 0, 1, 2, 3. 

U (x) = -\v{x), 
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Ui(x) = ±(3V(x) 2 -V"(x)), 

U 2 (x) = — (-10V{xf + 5V'(x) 2 + 10V(x)V"(x)-V {4) (x)) , 
U 3 (x) = y (35V(xf - 70V(x) V\xf - 70V(x) 2 V"(x) +2lV"{xf 
+28 V'{x) V (3) (x) + 14 V(x) \/ (4) (x) - V {6) (x)) . 

Now we can solve the first N relations of fl3~9|) in terms of the functions 

Uj{x) 

j'-i 

a N _j(x) = AN-kUj-k-iix) + A N _j, j = 1, . . . , N. (41) 

fc=0 

Inserting the above expressions into the last equation of ( |39D we get 

N 

A_! + J2 A N ^ k U N - k (x) = (42) 

fc=0 

(when deriving the above equation we take into account that by convention 
a_i(x) = 0). 

Equation (pE2"D is the necessary and sufficient condition for the Schrodinger 
equation (HD to be invariant with respect to the first-order operator (|37|) , 
whose coefficients are polynomials in e of the order N. But it is easy to see 
that given a condition ([42]), equation ([!]) admits operator of the form (|37|) , 
whose coefficients are polynomials in e of an arbitrary order N' > N. 

Indeed, the invariance conditions for the operator 

/ N' \ N 1 

\j=0 J 3=0 

have the form (|38|) with = N'. Coefficients aj(x) with j = N',N' — 
1, . . . , N' — N are given by formulae fl4"T|), where one should replace N by N', 
and the remaining coefficients by force of relation ([S) read as 

JV-l 

aj {x) = A jk U k (x) +Aj, j = 0, 1, . . . , N' - N - 1, 

k=0 

where A Q , . . . , A^'-n-i are arbitrary constants, Aj k are constants expressed 
via A N ,_ N , A N t. 
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Substituting these results into the last (j = 0) equation from (fjl|) yields 

N-l 

A_i + A C/ (z) + ^ A 0fe £4 + i(z) = 0, 

fc=0 

whence we conclude that, provided A-i = A = 0, A ok — 0, k = 0, 1, . . . , N— 
1, the invariance conditions are satisfied. 

Thus, if equation (^) is fulfilled with some N e N, then the corre- 
sponding Schrodinger equation admits arbitrary order Lie symmetries and, 
consequently, is exactly solvable. 

A general solution of the Schrodinger equation invariant under the oper- 
ator Q is given by (|33|). Substituting formula ( |33|) into equation flU), where 
V(x) is an arbitrary solution of (fE^), results in a (2iV + l)-th order alge- 
braic equation for e. Its solutions e = X±, . . ., e = X m are eigenvalues of 
the Schrodinger operator. Corresponding eigenf unctions are obtained if we 
insert e — Ai, ...,£ = A m into (p3|) . 

Summing up, we conclude that any solution V(x) of (^) gives rise to an 
exactly solvable Schrodinger equation. In what follows it will be established 
that the more strong assertion holds. Namely, if V(x) is a solution of (|42]) 
with some N, then the corresponding Schrodinger equation may have an 
arbitrary spectrum. 

Theorem 5 Let V = V(x) be a solution of ordinary differential equation 
( fjjjj with some fixed N 6 N. Then, the Schrodinger equation (fjj is exactly 
solvable and, moreover, the Schrodinger operator S = d 2 x — V(x) may have 
an arbitrarily prescribed spectrum. 

Proof. We will give the principal steps of the proof omitting technical details. 
As the potential V(x) satisfies equation (^2]), the corresponding Schrodinger 
equation admits a Lie symmetry of the form ([37]). Excluding from (|37|) the 
parameter e we recover symmetry operator Q of the order n = 2N + 1 
which commutes with the Schrodinger operator S. Next, we construct an 
operator Q\ = Q + f{s), where f(e) is an arbitrary smooth function. Evi- 
dently, Q\ commutes with S and, consequently, is a symmetry operator of the 
Schrodinger equation. This means that conditions of Theorem 1 are fulfilled 
and a fundamental system of solutions of ordinary differential equation 

Qiu(x) = 
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forms a basis of the invariant space V„ of the corresponding Schrodinger 
operator S. Representation of S in the space V n is given by the formulae 
fl2"5|), fl2"E|), where Ax = a + f\e) and parameters a, A 2 , ■ ■ ., A n , B%, . . ., B N 
are independent of f(e). Eigenvalues of the operator S are solutions of the 
characteristic equation for the matrix R having the entries (p5|), i.e. of the 
equation 

det \\Rjk — sSjk\\ = 0. 

It is not difficult to become convinced of that the above equation can be 
represented in the form 

Po(e)+ Pl (e)f(e)+ P2 (e)f(6) 2 = 0, 

where po, pi, p 2 are polynomials in e of the order not higher than n. 
Consequently, zeros of the function 

F(e)=po(e)+p 1 (e)f(s)+p 2 (s)f(e) 2 

are eigenvalues of the Schrodinger operator. 

As f(e) is arbitrary, the function F(e) may have an arbitrarily prescribed 
set of zeros and, thus, a spectrum of the initial Schrodinger equation ([[]) may 
be arbitrary. 

As an illustration, we will consider the simplest case when there exists 
such Ni that Un x (x) = 0. In such the coefficients of the symmetry 

operator fl3T|) with N = kN±, fceN are easily shown to be 

iH.r.c) = I J2 U i( x ) £Nl ~ j ~ l 

.7=-l 



-i 




where A , . . . , A^, B , . . . , are arbitrary constants. 

Inserting this result into ([IT]) and integrating we come to the following 
Ansatz for the function ip(x): 



ip{x) = Mb 1 Uj{x)e N iS- 1 



ax 
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where f(e) = (£f =0 B j£ ^ =0 A^)~\ 

Substitution of the expression (^) into the initial equation after some 
manipulations gives the following algebraic equation: 

N x -\ 

£ 2N 1+ i + j2 - /(e) 2 = 0, (44) 

j=0 

where Xj are integrals of the ordinary differential equation Un^x) = and, 
consequently, are constant on the set of its solutions. 

Thus, for an arbitrary Ni G N any solution V = V(x) of the ordinary 
differential equation Un 1 (x) = gives rise to an exactly solvable Schrodinger 
equation. Eigenvalues of the Schrodinger operator S are the roots of the 
algebraic equation (|44D and the eigenfunctions are of the form fl43|) , where e 



is an arbitrary solution of (|I3). 

It is instructive to consider in more detail the above formulae for the case 
N — N\ — 1. With this choice of iV the formula (|43D reads as 

Vi(x) = (2 e -V(x))3exp|-2/(e) (j ( 4 5) 

where f(e) = {Aye + Aq)(B\e + -Bo)" 1 and the function V(x) is a solution of 
the ordinary differential equation U±(x) = 0, i.e. 

- 3V(x) 2 + V"(x) = 0. (46) 

Inserting the Ansatz (|45| ) into (|T]) yields the following equality: 

-8e 3 +X + 8f(ef = 0, (47) 

where 

J = V(xf - l -V\xf 

is the first integral of equation ( |I6|) . Note that an alternative derivation of 
the formulae (|4^D-(|47D has been obtained in [P5[] . 



Now let us make an important remark. It is readily seen from formulae 
(H)~(I3) ^ na ^ function f(e) is not obliged to be a ratio of two first- 
order polynomials. It may be arbitrary. And what is more, eigenvalues of 
the Schrodinger operator are the zeros of the function F(e) = — 8e 3 + X + 
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8f(e) 2 . Choosing an arbitrary function f(e) in a proper way we can get the 
function F(e) having arbitrary prescribed set of zeros. This means that the 
Schrodinger equation (JTJ) with V(x) satisfying ( f4"T| ) may have an arbitrary 
spectrum. 



For example, if we choose / = ye 3 — X /8, then the function (|45|) is a 
solution of the Schrodinger equation under arbitrary e (the case of a con- 



tinuous spectrum). Next, if we choose / = ye 3 — X /8 + J2jLo Ajei, then 
a finite discrete spectrum is obtained (eigenvalues are roots of the N-th or- 



der polynomial). At last, choosing / = ye 3 — Iq/8 + sine yields an infinite 
discrete spectrum (eigenvalues are zeros of the sine). 
Similar results are obtained for N = N\ = 2 

tp(x) = (8s 2 -AeV(x) + 3V(x) 2 - V"(x))* 

dx \ 1 



xexp|-8/(e) I J 



8e 2 - AeV(x) + 3V(x) 2 - V"(x) 



where f(e) is an arbitrary function, V(x) is a solution of the ordinary differ- 
ential equation 

10V(x) 3 - 5V'(x) 2 - 10V(x)V"(x) + V i4) (x) = (48) 
and e is a solution of the equation 

-128e 5 - 2el x + J + 128 f(ef = 0. 
Here X , X\ are integrals of equation fl48|) of the form 
J = 5V(x) 4 -10V(x)V'(x) 2 - V"(x) 2 + 2V'(x)V {3) (x) 1 



I V{x) 19V{x) 5 J 2 51 V{x) 4 25V(x 



2 2 8V'(x) 2 AV'(x) 2 8V'(x) 2 

Z V"(x) 2 5V(x)V'(x) 2 V"(x) 4 



+ 



AV'(x) 2 AV'(x) 2 8V'(x) 2 ' 

Generically, if V(x) is a solution of the 2iV-th order ordinary differen- 
tial equation Un(x) = 0, then the corresponding Schrodinger operator S 
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may have an arbitrary spectrum. Eigenvalues of S are obtained by solv- 



ing the algebraic equation (|44D and its eigenf unctions by substituting the 



corresponding values for e into (H3T) . 



We will finish this section with one more puzzling property of the ex- 
actly integrable models obtained. Let us denote the total derivatives of the 
functions Uj(x) with respect to x as Wj(x) and consider an infinite set of 
evolution equations for a function u = u(t, x) 

^^ = F 3 [u(t,x)], j>l, (49) 

where the functions Fj are obtained from Wj(x) by formal replacement of 
V(x) with u(t, x). Now we see that equations obtained in this way form the 
famous integrable KdV hierarchy. Taking, for example, j = 2 yields the KdV 
equation 

du(t,x) 1 / , .du(t,x) d 3 u(t,x)\ 
—gf- = & { &U[t,X) ~dx dx^J ' 

Furthermore, differentiating the relations (^) we obtain the recurrence 
relations determining Wj(x) 



W j+1 (x) = YW j (x) = (-d 2 x -V(x)--V'(x)d^) W s {x), J = 0,1, 



with Wq(x) = 0. The operator Y above is nothing else but the well-known 
recurrence operator for the KdV hierarchy [j32]j . 

Next, if we formally replace V(x) by u(t,x) in Uj(x) determined by the 
recurrence relations (fl"0"D , then the densities of motion constants of the KdV 
equation are obtained, X being the recursion operator connecting these den- 
sities. 

Equations of the form (fE^) are known in the literature as the station- 
ary Lax-Novikov equations. Equations of the stationary KdV hierarchy are 
particular cases of equations ([EJ) with A\ — ■ ■ ■ — A^-i = 0. We have 
proved that any solution V = V(x) of the stationary Lax-Novikov hierarchy 



(42) yields an exactly solvable Schrodinger equation. The correspondence be- 



tween solutions of stationary KdV hierarchy and exactly solvable Schrodinger 



equations with reflection-less potentials is known (see, e.g. the paper |33j and 



references therein). Moreover, it has been established that the Schrodinger 
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operators with so chosen potentials may have an arbitrarily prescribed spec- 
trum. But the fact that solutions of the stationary Lax-Novikov hierarchy 
(fi2"D have the same property seems to be new. 



IV. Conclusion 



In view of numerous excellent papers and monographs (see e.g. [17] and 
the literature cited therein) devoted to developing algebraic methods for the 
investigation of spectral properties of equation (P, it is, of course, not enough 
to say that the problem of describing part of the spectrum of the Schrodinger 
equation is equivalent to computing its conditional symmetries in order to 
justify a necessity of introducing a new complicated structure. Our principal 
motivation for looking for a symmetry interpretation of the results obtained 
in this field is that it may open a possibility 

• to study the spectrum of two- and three-dimensional Schrodinger equa- 
tions, 

• to investigate 'spectral properties' of nonlinear Scrodinger equations, 
and 

• to study spectral properties of matrix differential operators (say, of the 
Dirac operator) 

by purely algebraic means. 

For instance, there are strong evidences that a necessary condition for a 
three-dimensional Schrodinger equation to be exactly solvable is an invari- 
ance with respect to a three-dimensional Lie algebra of high order symmetry 
operators. We can guess that one of the necessary conditions of 'quasi exact 
solvability' of the three-dimensional Schrodinger equations is a non-trivial 
conditional symmetry admitted. The simplest possibility to move in this di- 
rection is to combine the technique developed in the present paper with the 
method of separation of variables []34] . In our paper |3Sj we have classified po- 



tentials V(xi,X2) such that the corresponding two-dimensional Schrodinger 
equation 

ilpt + Ipxixi + 1px 3 x 3 = V(xi, x 2 )ip 

can be separated into three ordinary differential equations. One of these is a 
first-order equation and can always be integrated by quadratures. Two other 
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are exactly of the form (Jl|) and can be solved within the framework of the 
approach described in Sections 2 and 3. 

Furthermore, the method of conditional symmetries is applicable not only 
to linear partial differential differential equations but also to nonlinear ones 



21, 22]. Let, for example, the one-dimensional nonlinear Schrodinger equa- 



tion 

ip xx = (e + V(x) + V*, ^, €)) i> (50) 

be conditionally invariant with respect to an n-th order Lie-Backlund oper- 
ator 

Q = r ] {x,iP,iP',...,ipW)d i> + --- 



in the sense of |[22| . Then, using a technique similar to that developed in 
Section 2 we can construct an Ansatz for a function ip(x), which gives a solu- 
tion of (|50D , provided the energy parameter e satisfy some algebraic equation 
G{e) = 0. Solutions of this equation can be interpreted as eigenvalues of the 
nonlinear Schrodinger operator d x — V — F. 

An interesting example is a family of nonlinear Schrodinger equations 
suggested by Doebner & Goldin |36j. Taking the polar decomposition 

iP{t,x) =e r{t ' x)+ls(t ' x \ (51) 

and fixing the gauge V\ — — 1, v 2 = (this is always possible ||37|1 ) we can 
represent the Doebner- Goldin model in the following way: 

s t + 2[i 2 r xx + ii x s xx + 4(/i 2 + Hs)rl (52) 
+2(pi + fJLi)r x s x + fJLzS 2 x + HqV{x) = 0, 

where /xo, ■ ■ ■ ,^5 are model parameters. 

If we impose on the solutions of (|52"D an additional condition s x = 
(which picks out a subset of stationary solutions), then the system obtained 
is consistent if and only if the conditions 

r = r(x), s = —et, e = const G R (53) 

are fulfilled. 

With this choice of functions r and s system (|52|) reduces to a single 
equation 

- e + 2 M2 r" + 4(/i 2 + ^){r'f + fi V(x) = 0, (54) 
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which is either linear (/i2 + /is = 0) or can be linearized by the substitution 

^ {x)=exp {^h^ r{x) } (55) 

to become 

f" = ( 7T^ \ - 77-^ sV(x)) (p. (56) 

\4(ji 2 + (i 6 ) 4(/i 2 + /i 5 ) v >) Y K ' 

As established in Section 2, any equation of the form (|56| ) possesses a 
nontrivial high order conditional symmetry. Since the nonlinear equation 
(|5iD is equivalent to fl56|) , it possesses high order conditional symmetry as 
well. Thus, the initial Doebner-Goldin equation has a subset of solutions 
with non-trivial conditional symmetry which can be effectively applied to 
construct finite or even infinite (if the conditional symmetry can be reduced 
to a high order Lie symmetry) set of its exact solutions. 

We hope that the reasonings above are convincing enough to motivate a 
further study of high order conditional symmetries of linear and nonlinear 
Schrodinger equations in one, two and three dimensions. It may be also very 
interesting to study classical and conditional symmetries of the stationary 
Dirac equation in the presence of non- vanishing electro-magnetic field and to 
apply these to derive a spectrum of the Dirac operator. These problems are 
under investigation now and will be a topic of our future publications. 
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